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Abstract 

We use discrete holomorphic polynomials to prove that, given a 
refining sequence of critical maps of a Riemann surface, any holo- 
morphic function can be approximated by a converging sequence of 
discrete holomorphic functions. 



1 Introduction 

The notion of discrete Riemann surfaces has been denned in We are 
interested in discrete surfaces given by a cellular decomposition <) of dimen- 
sion two, where all faces are quadrilaterals. We suppose ^> bipartite and it 
defines (away from the boundary) two dual cellular decompositions T and 
T*, edges are dual to edges Ti, faces are dual to vertices T and vice- 
versa. Their union is denoted the double A = T U T*. A discrete conformal 
structure on A is a real positive function p on the unoriented edges satisfying 
p(e*) = l/p(e). It defines a genuine Riemann surface structure on the dis- 
crete surface: Choose a length 5 and realize each quadrilateral by a lozenge 
whose diagonals have a length ratio given by p. Gluing them together pro- 
vides a flat riemannian metric with conic singularities at the vertices, hence a 
conformal structure ||. This data leads to a straightforward discrete version 
of the Cauchy- Riemann equation. A function on the vertices of <0> is discrete 
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holomorphic iff for every quadrilateral (x,y,x',y') G ^2, 

f(y')-f(y) = i P (x,x') (f( x ')-f( x )) 



(l.l) 



9 y> 

>£(y,y') 



> 



X 





Figure 1: The discrete Cauchy-Riemann equation. 



Given a fixed flat riemannian metric on a Riemann surface, with a finite 
number of conic singularities, we define a discrete conformal structure as 
critical if the flat riemannian metric it gives rise to is isometric to this fixed 
one. Let's stress that this notion is useful mostly when considering a sequence 
of discrete conformal structures, all adapted to the same fixed flat metric. 
A refining sequence is such a sequence of critical maps where the common 
side lengths 5k of the diamonds in the map ()k, form a sequence converging 
to zero, and where the lozenge angles are bounded away from zero (the faces 
don't collapse). 

In ||1| we proved that 1) any Riemann surface accepts a refining sequence 
of critical maps and 2) a converging sequence of discrete holomorphic func- 
tions on a refining sequence of critical maps converges to a continuous holo- 
morphic function. It is the main purpose of this article to prove the converse 
result: 

Theorem 1.1 Given a refining sequence of critical maps of a compact flat 
simply connected surface U, any holomorphic function on U can be approxi- 
mated by a converging sequence of discrete holomorphic functions. 
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The proof is based on series expansion. In Sec. ^| we define discrete inte- 
gration at criticality and show that the sequence of primitives of converging 
discrete holomorphic functions converge to the continuous primitive. We 
define discrete polynomials and use them to prove the main theorem. 

In Sec. H, we define discrete derivation and present its differences with 
the continuous one. It allows none the less for the coefficients of a discrete 
polynomial to be computed by successive derivation. In the annex, we discuss 
different aspects of the discrete theory which differ or don't appear in the 
continuous case, in Sec. |A] another notion of derivation, in Sec. [B] the analog 
of the Leibnitz formulae, in Sec. |C| a motivation of the derivation formula 
through a study of the minimal polynomial and in Sec. [D] a discussion of the 
problems arising when considering discrete series, with the example of the 
discrete exponential. 



2 Integration at criticality 

Away from a conic singularity, a critical map locally forms a flat partition of 
the plane by lozenges. The relevance of this kind of maps in the context of 
discrete holomorphy was first pointed out and put to use by Duflin 0. He 
defined the discrete analogues of the integer powers of Z and the derivation 
formula that we are going to give. 

The crucial point about criticality is the following remark: Given an iso- 
metric local map Z : U H A — > C, where the image of the quadrilaterals are 
lozenges in C, any holomorphic function / G Q(A) gives rise to an holomor- 
phic 1-form fdZ defined by the formula, 

fdZ := (Z(y) - Z{x)) , (2.1) 

where (x, y) G is an edge of a lozenge. It is averaged into a holomorphic 
1-form on A: 

fdz= m±m+m±m (z(x>) . z(x)) . (2 . 2) 

(x,x>) 4 

Alternatively, it provides a way to integrate a function by taking the 
primitive of fdZ. 
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Theorem 2.1 Given a sequence of discrete holomorphic functions (f k ) on 
a refining sequence of critical maps, converging to a holomorphic function f , 
the sequence of primitives (J fkdZ) converges to J f{z)dz. Moreover, in the 
compact case, if the convergence of the functions is of order 0($fy, it stays 
this way for the primitives. 

Proof \2. 1\ . We recall from [jlj that we extend a discrete holomorphic 
function fd from the vertices <Oo linearly to the edges <0i and harmonicly to 
the faces <^2 to obtain a continuous piecewise harmonic function fd of the 
surface. We proved that, given a refining sequence (<0 fc ) of critical maps, the 
point-wise convergence of a sequence of discrete holomorphic functions (fk), 
restricted to the sequence of vertices (<)($), implies a uniform limit of the 
continuous extensions (fk) to a genuine (continuous) holomorphic function 
/• 

Let U be the flat simply connected patch under consideration. We sup- 
pose that we are given a sequence of flat vertices Ok G <$k where the face 
containing the fixed flat origin O G U is adjacent to Ok- Let Fk the extension 
of the primitive f Q fk dZ to the whole surface. We want to prove that the 
following sequence tends to zero 



{F k {x) - F k {0)) 



f(z) dz 



o 



(2.3) 



ken 



Let us fix a point x, and for each integer k consider a vertex Xk G ^>o on 
the boundary of the face of <0>2 containing x. 

We decompose the difference ( ^3|) into three parts, inside the face con- 
taining the origin O and its neighbor Ok, similarly for x and Xk, and purely 
along the edges of the graph <) fc itself. 

\(F k (x) - F k (0))- f f(z)dz\ = 

Jo 

(F k (x) - Fk(x k )) + [ " fkdZ + (F k (0)-Fk(Ok))- [ f(z)dz 
Jo k Jo 



< 



F k {x)-Fk{x k ) - \ f{z)dz 



Xk PXk 

fkdZ - / f(z) dz 
o k Jo k 



+ 



F k {O k ) - F k (0) 



o k 



f(z) dz 



o 



(2.4) 



On the face of <0 containing x, the primitive x i— > f(z) dz is a holomorphic, 
hence harmonic function as well as x i— > F k (x). By the maximum principle, 
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the harmonic function x \— > F k (x) — F k (x k ) — f(z) dz reaches its maximum 
on that face, along its boundary. The difference of the discrete primitive along 
the edges of <0> is equal to 



Fjfe((l - X)x + \y) - F k {x) = \{y - x) 



fk{x) + f k (y) 



(2.5) 



while because / is differentiate with a bounded derivative on U, 



(1-X)x+Xy 



f{z)dz = X(y-x) 



f(x) + f(y) 



+ (y-x) 



,\*f'(x) + (i-\ynv) 



X(y-x)M+M + 0{ 5l) 



+ o(5l) 
(2.6) 



so that 



\F k (x) - F k {x k ) 



f{z)dz\=0{5l) 



(2.7) 



Similarly for the term around the origin. 

By definition of f k , the 1-form f k {z)dz along edges of the graph <0> is 
equal to the discrete form f k dZ so that f Q k f k dZ = J** f k (z) dz on a path 
along ^> edges. Therefore the difference 



Xk r^k 

f k dZ- / f{z)dz 
o k Jo k 



< 



(fk(z) - f(z] 



dz 



(2.8) 



is of the same order as the difference \f k (z) — f(z)\ times the length £(7^) of 
a path on <>& from O k to x k . We showed in 0] that provided that the lozenge 
angles are in the interval (77, 27r — 77) with 77 > (the faces don't collapse), this 
length can be bounded as £(%) < — Ok\- Since we are interested in the 

compact case, this length is bounded uniformly and the difference ( |2.8| ) is of 
the same order as the point-wise difference. We conclude that the sequence 
of discrete primitives converges to the continuous primitive and if the limit 
for the functions was of order 0(5 2 ), it remains of that order. <0 

Following Duflin |4], [|, we define by inductive integration the discrete 
analogues of the integer power monomials z k , that we denote Z 



k:. 



Z :0: 
Z :k: 



Z :k ~ 1: dZ, 



(2.9) 
(2.10) 



o 
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where O is a fixed flat origin. The choice of this origin is relevant and we 
present in the Appendix the formulae for the change of basis of discrete poly- 
nomials from one base point to another. Although polynomials are translated 
into polynomials of same degree, the formulae are slightly more involved than 
the Leibnitz rule (Z - a) n = J2l= ©^ fc a n ~ fc . When sequences of maps are 
involved, we suppose that we have chosen a flat origin for each map, forming 
a converging sequence to a flat point. For example we suppose that a fixed 
flat point O on the surface is a vertex of each map. 

The discrete polynomials of degree less than three agree point-wise with 
their continuous counterpart, Z :2 '(x) = Z(x) 2 . 

A simple induction then gives the following 



Corollary 2.2 The discrete polynomials converge to the continuous ones, 
the limit is of order 0(51). 

Which implies the main theorem: Proof [Jj| . On the simply connected 
compact set U, a holomorphic function / can be written, in a local map z as 
a series, 



f(z) = J2*kZ k . (2.11) 



feeN 



Therefore, by a diagonal procedure, there exists an increasing integer se- 
quence (N(n)) nm such that the sequence of discrete holomorphic polynomi- 
als converge to the continuous series. 




(2.12) 







3 Derivation at criticality 

We are interested in getting the coefficients of a polynomial by the same 
process as in the continuous case: successive derivation. We are going to see 
that this operation, while possible, is more difficult than in the continuous. 

Let e be the biconstant e(T) = +1, s(T*) = —1. The main issue is to 
separate this "mode" from the constant mode in order to define properly the 
discrete analog of the value of a function at the origin. 
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On a finite map, the dimension of the space of discrete holomorphic func- 
tions is equal to half of the number of boundary points plus one, 

n(u) ~ c |ai/|/2+1 , (3.1) 

for example the values on the boundary points Tq H dU on the graph T* and 
in one point O G T. The dimension of the space of discrete polynomials can 
not be greater than this number. This implies the existence of a minimal 
polynomial 

n 

P z = J2^Z- k: = 0. (3.2) 

k=l 

For a holomorphic function /, the equality / dZ = is equivalent to / = Xe 
for some A G C, we can uniquely normalize Pz so that Yuk=i ^ ak Z'- k ~ h = e, 
that is ai = 1. 

We can therefore take as a basis of the discrete polynomials the set 

[l.Z.Z^ Z» \ :). (3.3) 

We conjectured in JIJ that when U is convex [[J, ( |3.3| ) is a basis of the 
whole space of discrete holomorphic functions. In any case, one can supple- 
ment them by an orthogonal complement and we will denote v G & 9U ^ 2+1 
the vector encoding the linear form yielding the coordinate along the constant 
function 1 in this basis: 



n-2 



v ■ a k z :k - 



a , 



k=0 



v-e = 0, (3.4) 
v ■ f = when VJfe, Z :k - ■ f = 0. 

Notice that since the change of base point for polynomials is triangular in 
the basis of monomials, the vector v does not depend on the map Z or on the 
base point but only on the discrete conformal structure. An explicit value 
for this vector would be desirable, especially for numerical purposes. 
Following Duffin j|, [| (see Sec.|C| for a motivation), we introduce the 

Definition 3.1 For a holomorphic function f , define on a flat simply con- 
nected map U the holomorphic functions f \ the dual of f , and f, the deriva- 
tive of f , by the following formulae: 

/+(*):= £(*)/(*), (3-5) 
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where f denotes the complex conjugate, e = ±1 is the biconstant, and 

f'(z)--=^(£pdZ^ +\e, (3.6) 

with A determined by 

A = v ■ /'. (3.7) 

If, in the context of sequences, the function / can be Taylor expanded around 
the origin, f(x) = f(0) + fix + 0(x 2 ), then A = /i and the identity [TT] is 
equivalent in the 0(x 2 ) sense to 

E^ g *)( z^"^ - A ) =0 ' (3 - 8) 

where (0,yk,Xk,yk+i) £ ^2 are the quadrilaterals adjacent to the origin. 
Therefore this condition can be used for numerical purposes in replacement 
of the exact derivation formula |3.7[ It states that the derivative at the origin 
is the mean value of the nearby face derivatives (see Sec. 0). For example, 
in the rectangular lattice <0> = <5(Ze i61 + Ze - * 6 '), with horizontal parameter 
p = tan# and vertical parameter its inverse, it provides the good choice for 
the first three and all the even degrees, A = 0, Ai = 1 and X 2 k = 0, but 
yields A& = k\ (0 cos(/c — 1)0 = 0(5 k ~ 1 ) for k odd and fixed, instead of 
for k > 2. 

Proposition 3.2 The derivative f fulfills 

do/ = f'dZ. (3.9) 

The discrete monomials verify 

vjfe<n-i, (z :k: y = kz :k ~ l: . (3.10) 

Eq.fl3.9|) was proved in [TJ and Eq.( |3.1U| ) follows from 



fdZ = <£> f = Xe. (3.11) 
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A Face derivation 



Another derivation operator can be denned. By definition of holomorphy, for 
each face of the graph <}, the 1-form d\f is proportional to dZ along its pair 
of dual diagonals. This defines a linear operator 

^■M(Anu)^c\Onu). (A.i) 

Given the local flat map Z, the space C 2 {() D U) can be seen as the space 
(7(1.0) ^ n u) of type (1, 0) forms on A, the -M-eigenvectors of the Hodge star 
*, that is to say the forms a such that, for dual edges (y, y') = (x, x')*, 

a = ip(x,x') / a. (A.2) 
(y,y') J(x,x') 

It defines 

a z EC 2 (On U) for a = a z dZ (A.3) 

We can say that the face function az is holomorphic whenever the corre- 
sponding form a is closed, that is to say by imposing the Morera theorem. 
This idea was developed by Kenyon |J. Notice that the 1-form a is then co- 
closed as well. If the quadrilateral graph <0, along with its discrete conformal 
structure was isomorphic to its dual graph <)*, a choice of isomorphy could 
define a derivation endomorphism ^ : fi(A D U) — > f2(A fl U). But for this 
to happen, the vertices have to have degree four, so although important, the 
only case is Z 2 with the choice of a translation by (n + ^, m + |), n, m G Z. 

It is immediate that given a 0(<5 2 )-converging sequence of discrete holo- 
morphic functions, this face derivation yields a sequence of functions, ex- 
tended on each face as piece-wise constant functions, which converges to the 
continuous derivative. 



B Change of basis 

We are now going to consider the change of coordinate. Let Z a critical 
map. If C is another critical map with Q = a(Z — b) on their common 
definition set, the change of map for Z :k: is not as simple as the Leibnitz rule 
(a (z — b)) k = a k X^Lo ^^'^{^y but is a deformation of it. The problem 
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is that pointwise product is not respected, Z :k+t (z) ^ Z :k '\z) x Z :t (z). In 
particular the first is holomorphic and not the second. For each partition of 
m = k + £ into a sum of integers there is a corresponding monomial of degree 
m. An inference shows that the result is still a polynomial in Z: 

Proposition B.l The powers ( [k: of the translated critical map ( = a(Z — b) 

are given by 



a k (-l) j Z tk - j: B j (b) (B.l) 



= u 

j=0 



where (b) corresponding to IP is a sum over all the degree j monomials in 
b, defined recursively by B° = 1 and 

B k (b) := (-l) k+j+1 Z :k - j: (b)B j (b) (B.2) 



Proof \B7l\ . The calculation is easier to read using Young diagrams. Denote 
the point-wise product of monomials 

[Z^\z)) i2 . . . (Z-- k -\z)) ln (B.3) 

with k\ > hi > . . . > k n as a Young diagram Y, coding columnwise the 
partition of the integer given by the total degree k = YTj=i kj x 4» the 
sum of t = YTj=i integers. For example the following monomial of degree 
15 = 3x2 + 2x4 + 1 is noted 



(Z^(z)Y (Z^z)y Z-iz) =: ^ffl (B.4) 



Then, B^(b) = X]y c (^0^(^) where the sum is over all Young diagrams 
of total degree j, c(Y) is an integer coefficient that we are going to define 
and Y(b) is the pointwise product of the monomials coded by Y at b. The 
coefficient of the Young diagram Y above is given by the multinomials 



r(Y) = (-l) k+i - - 



o 
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B° 
B 1 
B 2 
B 3 

B 4 
B 5 

B 6 



+ 1 
+ □ 

-] + 2m 

-8P + 6g-36^ + 24 Em 



+ 



- 10 



- 20 |3 + 60 ^ + 90 gp - 240 + 120 

1 + 12 FT + 30 =3 - 90 1° + 20 g - 360 gP + 480 J 
90 + 1080 - 1800 + 720 



i 



Table 1: The first analogs of z k needed in a change of basis. 



For example, the first Young diagrams have the coefficients c( i i a i i ) 



9,! T?.! 



1!-1! re! 



77.!, C( 



ID ) 



) = ("I) 

(n+l)! 



n+1 



2! l!(n-l)! 



= (-ir^ifTi = (-ir +l 2(n+i) 

. (rH-l)! n. ^ g rg ^ ^ ermg are listed 



and c(EJ 

explicitly in Table [l|. 

It is to be noted that the formula doesn't involve the shape of the graph, 
the integer coefficients for each partition are universal constants and add up 
to 1 in each degree. As a consequence, since in the context of a refining 
sequence of critical maps, 



Z :k: (z)Z :e: (z) = Z- k+t (z) + 0(S 2 



(B.6) 



with k,£,z fixed, the usual Leibnitz rule is recovered in 0(5 2 ). Let's stress 
again that these functions B k are discrete functions on the graph A which 
are not discrete holomorphic. 



C Minimal polynomial 

We give a motivation to Duffin's derivation formula Eq. ( |3.6|) by studying 
the minimal polynomial, its behavior regarding duality t and its derivatives. 
It shows that derivation in Z acts as integration in Z\ 
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Consider the minimal polynomial 

n 

p z = J2^z- k - = o. (C.l) 



k=l 



with d\ = 1, or equivalently J2k=i ^ a k Z' k 1: — £ — ^ ■ This implies that 

n 4 
J2Hk-l)a k Z-- k ~ 2 -- = -& + a 2 2e, (C.2) 



k=2 



5 2 

2 7 :2:t 



n / 4 \ 

= $> n+1 _J-J 

fc=0 ^ ' 



4 V(n + l-fc)! z:fc t_ 



fc! 

Since Pz(Zy = P Z (Z^), a symmetry among the coefficients follows: 

afc = ^^ — ^rUJ ' (a3) 

leading to 

u-2fc + l ra-2(ra+l-fc) + l 

fc! l afc l(^) * = (n + 1 - k)\\a n+1 _ k \ (J^j 4 , (C.4) 

n-l 

in particular |a n | = (^-) 2 and when n is odd, arg (on+i) = |arg (a n ). 

D Series 

We present here some problems linked to discrete series. They can converge 
only when their coefficients decrease at least exponentially fast. And expo- 
nentially fast is enough since one can define a discrete exponential. 

The 0(5 2 ) convergence of discrete polynomials to continuous ones is in 
fact slower and slower as the degree of the polynomial grows. Consider for 



12 



k 


3 


4 


5 


6 


7 


Z k (x) 


™3_i_ x 


n z 


5 _|_ 5x a _j_ 3x 
X "r n 2 I" 2n 4 


6 1 10x 4 | 23x 2 


7 1 35x 5 | 49z 3 | 45ic 


n 2 n 4 


n 2 n 4 4«6 



Table 2: The first powers on the interval [0, l]/n for x 



'n. 



example ^> containing the chain {0, . ,1}, the first Z :k '\x) with x = - 

are listed in Table |^. And in general, a scaling argument shows that 

\Z'- k '-(x) -x k \ < X k \x\ k - 2 5 2 . (D.l) 

But the constant A& is growing exponentially fast with k. It is for example 
not true that, at a point close enough to the origin, Z :k: will tend to zero 
with increasing k. On the contrary, if x = b e %e is a neighbor of the origin 
with (0,x) G ()i and k > 1, then 

^(*) = — (D.2) 

in fact diverges with fc. If y is a next neighbor of the origin, with the rhombi 
(O, x, y, x') G <^2 having a half angle 9 at the origin, 

has the same diverging behavior, and so has every point at a given finite dis- 
tance of the origin. It is only in the scaling limit with the proper balance given 
by criticality that one recovers the usual behavior \x\ < 1 ==>- \x k \ > 0. 

Therefore, for its general term to converge to zero, the coefficients of a 
discrete series must decrease at least exponentially fast. The discrete expo- 
nential Exp (:A: Z), for |A| ^ 2/5, is defined by 

Exp(:A:0) = 1 (D.4) 
dExp(:A:Z) = XExp (:X: Z) dZ (D.5) 

This discrete holomorphic function was first defined in [Q and put to a very 
interesting use in |J. The discrete exponential is a rational fraction in A at 
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every point, 



k 2 C 

where (9 k) are the angles defining (5e l9k ), the set of (Z-images of) <C>-edges 
between x = ^<5e l6lfc and the origin. Expanding exp(log(Exp (:A: x))) in 
0(5 2 ), for |A| < 2/5 and x fixed in a refining sequence of critical maps, we 
get that 

Exp(:A:x) = exp(Ax) + 0(5 2 ). (D.7) 



For |A| < 2/5, the series ^2^=0 A 5' is equal to (p.6|) whenever the former is 



defined, as its exterior derivative fulfills the right equation (p.5| ) and its value 
at the origin is 1. The great difference with the continuous case is that the 
series is absolutely convergent only for bounded parameters, |A| < |. This 
suggests that asking for a product such that 

Exp (:A: Z) ■ Exp Z) = Exp (:A + p. Z), (D.8) 

or equivalently Z :k - ■ Z :t = Z'- k+t may not be the right choice. The symmetry 
f is interesting as well, 

Exp(:A:Z) t = Exp(:i:Z), (D.9) 
A 

in particular, Exp (:oo: Z) = e. 

The general change of basis of a given series however possible in theory is 
nevertheless complicated and the information on the convergence of the new 
series is difficult to obtain. The exponential remains a particularly simple 
case, if C = a (Z ~ b)'- 

fc=0 ' fc=0 
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